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■ The derivation of the viscosity formula of accretion disks appearing in a number of 

' textbooks is based on a mean free path theory of gas particles. However, this procedure, 

when foUowed precisely, leads to the incorrect conclusion that the angular momentum of an 
accretion disk flows from the outer part to the inner part. 
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Introduction The standard model of accretion disks, also called the a disk model, 
proposed by Shakura and Sunyaev^^ states that the angular momentum is trans- 
ported from the inner part of the disk to the outer part, due to the action of some 
I kind of viscosity. The Reynolds number of the flow in accretion disks is estimated to 

^!t| ' be as high as 10^^ , which means, according to the accepted hydro dynamical wisdom, 

that the flow is turbulent (Spruit^)). It has been assumed that in calculating the 
viscous stress due to the turbulence one can employ the usual molecular viscosity 
formula. 

. There are many textbooks and reviews about the formation of accretion disks 

Oh| and the mechanism involved in the transport of angular momentum (e.g., Pringle, ^\ 

Q ■ Prank, King and Raine,^^ Spruit,^-* Hartmann, and Kato, Fukue and Mineshige, 

^^). Some of these textbooks, however, give elementary explanations about the 
^ I phenomenon of angular momentum transport based on physical arguments, i.e. a 

mean free path theory, rather than with mathematical rigor. This is perhaps because 
the authors of such textbooks intend to make the explanation easier for beginners 
^ ' to understand. We find that these explanations in fact do not lead to the results 

d • intended by the authors but to the incorrect conclusion that the angular momentum 

in an accretion disk is transported inwards. To the contrary, the angular momentum 
must be transported outwards, or accretion will not occur. There are, however, 
some other textbooks and reviews that give rigorous mathematical derivations of 
formulas describing the angular momentum transport, starting from the formula for 
the molecular viscosity (e.g.. Spruit, Kato et al. ^•*). 

The above mentioned "elementary explanations" based on the mean free path 
theory, however, seemed very plausible to us. We therefore began to doubt if the 
"molecular viscosity formula" itself is applicable to gas in an accretion disk rotating 
with Keplerian motion in the gravitational field of a central compact star. Eventu- 
ally, we found that these "elementary explanations" lead to the incorrect conclusion 
and that the "molecular viscosity formula" is basically correct. 
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Is mean free path theory applicable to the present problem? In the standard disk 
model, the equation for angular momentum conservation is 

where 

R = distance from the center of the central star, U = surface density, 
f2 = V(f,/R, and ly = kinematic viscosity. 

How can this equation be derived? It can be obtained using the formula of 
molecular viscosity for incompressible gas (e.g., Spruit 

/ dvi dvk\ 

where is the stress tensor, n is the viscosity coefficient, and Vi and are averages 
of the molecular velocity in the i- and A;-direction, respectively. The viscous force is 

Soifc ^ / _^ Ai;- (3) 

dxk dxk \dxk dxiJ 
A problem emerged when we were learning with a rather elementary textbook, which 
gives explanations for beginners. Since the situation is best described in the text- 
book, we make a rather lengthy quotation from Accretion Processes in Star Forma- 
tion, by Hartmann,^) p. 30 (the latter part of §5.1): 

Following Frank ct al. (1992), we can calculate the magnitude of the 
angular momentum transfer in terms of a kinematic viscosity. The basic 
picture is one in which turbulent elements of the gas moving at a typical 
random velocity w travel a mean free path A before mixing with other ma- 
terial. Thus, the net torques or angular momentum transfer at cylindrical 
radius R (Figure 5.1) will be produced (schematically) by the differing an- 
gular momenta of two streams of material; one from material originating 
at i? — A/2 and moving outward across R to mix with annular material 
centered at R+\/2\ and the other starting at R + X/2 and moving inward 
across R to mix with the inner annulus at -R — A/2. 

In this kinematic viscosity model, no net angular momentum is trans- 
ported unless there is shearing orbital motion, dfi/dR = i7' 7^ 0. The net 
angular momentum fluxes can be calculated in the following way. Material 
originating at i? — A/2 has an angular momentum of 

= {R- \/2fQ{R - A/2) = {R- \/2f[Q{R) - {X/2){dn/dR)], (5.9) 

where we have approximated the difference in the angular velocities in 
terms of the first derivative with respect to R. A similar expression with 
the negative values changed to positive applies to the material at A/2. 
The inner material diffuses outwards at velocity w and the outer material 
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diffuses inward at w across R. (The net inward motion of material in the 
accretion disk is assumed to be small in comparison with the turbulent 

velocity w). For simplicity wc integrate or average the disk structure in 
the direction perpendicular to the disk plane; then the net outward transfer 
of angular momentum across R per unit length for a disk with mass density 
per unit area E is 

Uw[{R-X/2f{-X/2)dn/dR-{R+\/2f{\/2)dO/dR] = -EwXR^dQ/dR, 

(5.10) 

where we have assumed that A is a short distance compared with the 
scale over which Q varies significantly. With this result, the total angular 
momentum flux outward across R, i.e. the torque of the inner annulus on 
the outer annulus, can be written as 

g = -27rRUi'vR^dn/dR, (5.11) 

where the viscosity is 

z/„ = \w. (5.12) 

Note that a negative gradient of angular velocity (i.e.i? decreasing out- 
ward) leads to a positive outward flux of angular momentum, as predicted 
by the qualitative picture of friction between neighboring annuli discussed 
above. 

Hereafter, we explain why the above derivation of (5.11) is not correct. In Eq. 
(5.9), {R — A/2) in the original textbook should be {R — A/2)^, and is so corrected 
in the present text. From (5.9) and a similar expression for Jout, we have 

Jin - Jout = f2{R)[{R - \/2f -{R + X/2f] 

- {X/2)df2/dR[{R - X/2f + {R + X/2f] 

= f2{R){-2RX) - {X/2)dQ/dR{2R^) = -Xd{R^Q)/dR. (4) 
Then, this gives, not Eq. (5.11), but 

g = -2-KRUu^d{R^Q)/dR. (5) 

The difference between Eqs. (5.11) and (5) is that the factor R? is operated on by 
the d/dR operator in Eq. (5) and not operated on by it in Eq. (5.11). If Eq. (5.11) 
is correct, the angular momentum flows in such a direction as to realize dQjdR = 0, 
that is, to produce rigid body rotation. In that case, the angular momentum flows 
outwards, which is what the standard theory predicts. 

Attention should be paid here to the following. That is, with a Keplerian disk, 
we have f2 oc = Rf2 oc and L = Rv^f, oc -R^/^, where is the 

circumphcrential velocity and L is the angular momentum. An accretion disk has 
the characteristic features that Q and decrease with increasing R and that L, in 
constrast, increases with increasing R. (See Fig. 1 to better understand the above 
description.) 
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If, in contrast, Eq. (5) holds, the angular momentum flows in the direction to 
realize d{R^f2)/dR = 0, that is, to a constant angular momentum. In this case the 
angular momentum flows inward. This means that even an originally disk-shaped 
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Fig. 1. Angular momentum transport between adjacent annuli in a thin accretion disk. The angular 
velocity Q and circumferential velocity — RQ decrease with increasing distance R. The 
angular momentum B?Q, however, increases with increasing R. The big wide arrows indicate 
which way the angular momentum is transported, inwards or outwards. 

gas will end up as a ring. We conclude here that the correct manipulation of Hart- 
mann's procedure leads to the conclusion that he did not intend. 

There is another famous textbook entitled Accretion Power in Astrophysics writ- 
ten by Frank, King and Raine.^^ There, the authors give a similar elementary ar- 
gument, with a similar confusion, while referring to a similar figure, as follows (on 
pages 58 and 59 of the second edition)*) : 

... As these elements of fluid are exchanged, they carry slightly differ- 
ent amounts of angular momentum: elements such as A will on average 
carry an angular momentum corresponding to the location i? — A/2, while 
elements such as B will be representative of a radial location R + A/2... 

If the authors had proceeded in this way correctly, they should have obtained 

*' Interestingly, this description has been revised and differs from the corresponding one in the 
first edition (1985). There, the "correct" conclusion is derived from the incorrect assumption that 
a gas element in a disk would, in chaotic motion, move radially a distance of order A away, while 
carrying its momentum (not its angular momentum). Since the rotating gas element should carry 
its angular momentum on that occasion, this assumption is wrong. This may have been the reason 
for the revision in the second edition. 
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the same result as from Hartmann's "correct" (i.e. corrected in the present paper) 

manipulation. They, instead, continue as follows: 

... As seen by an observer corotating with the fluid at P [at R] with an- 
gular velocity Q{R), the fluid at i? — A/2 will appear to move with velocity 
{R — \/2)Q{R — A/2) + Q{R)\/2. Thus the average angular momentum 
flux per unit arc length through R = const in the outward direction is 

pvH{R - \/2)[{R - \/2)Q{R - A/2) + Q{R)\/2] (a) 

[where v is the speed of gas elements in a chaotic motion, and the formula 
numbers (a) and (b) are temporarily put by the present authors]. An 
analogous expression, changing the sign of A, gives the average inward 
angular momentum flux per unit arc length. The torque exerted on the 
outer ring by the inner ring is given by the net outward angular momentum 
flux. Since the mass flux due to chaotic motions is the same in both 
directions, one obtains to first order in A the torque per unit arc length as 

-pvHXR^f}', (6) 

where i7' = dQjdr [and tl is the half thickness of the disk]... 

Although the meaning of (a) is not clear and is discussed below, let us accept (a) 
for the moment. We denote the formula (a) by Lj„ and make a linear approximation, 
leading to: 

Lin = pvH{R - \/2)[{R - \/2){f2{R) - {\/2)Q{R)') + n{R)\/2] 

= pvH{R - X/2)[Rf2 - R{X/2)f2']. (6) 

Similarly, 

Lout = pvH{R + X/2)[Rf2 + R{X/2)f2']. (7) 
Then, the torque should not be, —pvHXR^f?', but rather 

Lin - Lout = pvH[{R - X/2){RQ - R{X/2)f2') - (i? + X/2){RQ + R{X/2)n')\ 

= -XpvHR{f2 + Rf2') = -uERd{Rf2) /dR. (8) 

If this is true, the angular momentum flows in a direction to realize d{RQ)/dR = 0, 
which implies a constant azimuthal velocity and is unrealistic in a rotating gas, 
although in this case the angular momentum flows outward. 

Now we return to the problem of (a). When seen from P, the fluid at R — X/2 
should appear to move with velocity {R-X/2)Q{R-X/2) -RQ{R) + Q{R)X/2, which 
completes the transformation from the inertial frame to the corotational frame with 
the origin at P. Eq. (a) then should be 

pvH{R - X/2)[{R - X/2)Q{R - X/2) - RQ{R) + Q{R)X/2], (a)' 

so that 

Lin = pvH[{R - X/2){-R{X/2)Q') = -pvHR^{X/2)f2' , (6)' 
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L^n - L,ut = -pvHB?\Q'. (8)' 

This is the correct result. We then, however, notice that the procedure to ob- 
tain (8)' is simply that to obtain the well-known viscosity formula (2) , thus being no 
longer based on the mean free path theory. 

Discussion Having pointed out the confusion in the derivations given in the 
above cited textbooks, we now clarify how naturally the wrong conclusion comes 
from application of the mean free path theory (for discussion of this theory sec, for 
example, Vincenti and Kruger ^) Chap. 1, §4.) to a rotating gas in an accretion disk. 

A gas element rotating in a Keplerian orbit possesses an angular momentum, 
linear momentum and angular velocity. The element is then assumed to travel across 
R over a distance of A, the mean free path. If we denote the quantity to be trans- 
ported then by the element as Q, which is the angular momentum, linear momentum 
or angular velocity, then, by linear approximation, 

Q{R + A) - Q{R) = XdQ{R) /dR. (9) 

Thus, if we want to obtain the correct result, i.e. to introduce the angular velocity 
alone in d/dR, we have to make the impossible assumption that Q is the angular 
velocity. This argument clearly shows the failure of the application of the mean free 
path theory. 

Do we then have to use the molecular viscosity formula to consider the present 
problem? This is the next question. In a strict sense, the molecular viscosity formula 
(2) is derived for a non-rotating gas. A viscosity formula for a rotating gas should 
be derived using the Boltzmann equation with the Coriolis force taken into account. 
For a rigidly rotating gas, this was formulated by Chapman and Cowling, who 
considered the effect of the Lorentz force rather than the Coriolis force, with the 
result that the viscosity coefficient along the direction parallel to the rotation axis 
is the same as that in a non-rotating gas, while that perpendicular to the axis is 
suppressed in a manner that depends on the Knudsen number. The viscosity formula 
for a Keplerian rotating gas is much more complicated, and an exact formula cannot 
be obtained. Approximate formulas obtained by a few authors are summarized in 
the textbook of Pridman and Gorkavyi. 

The authors thank Professor S. Mineshige for his valuable comments and advice. 
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